Nanosecond Dynamics of Single-Molecule 
Fluorescence Resonance Energy Transfer 

Q : G.O. Ariunboldi'2'5, G.S. Agarwal^'^'^ Z. Wang^ 

C\i ; M.O. and H. Walther^ 

^ ■ ^ Max-Planck-Institut fiir Quantenoptik 

^ ■ D-85748 Garching, Germany 

^Institute for Quantum Studies, 
Department of Physics 
^ ' Texas A&M University, College Station, 

> : TX 77843-4242, USA 

^ ! ^ Physical Research Laboratory, 

Navrangpura, Ahmedabad-380 009, India 
^ Department of Chemistry, Princeton University 
^ ; Princeton, NJ 08544, USA 

O I ^ Theoretical Physics Laboratory, 

\ National University of Mongolia 

210646, Ulaanbaatar, Mongolia 



0^ 



O 



§ ■ February 1, 2008 

Abstract 

Motivated by recent experiments on photon statistics from individ- 
' ual dye pairs planted on biomolecules and coupled by fluorescence reso- 

5— ( ' nance energy transfer (FRET), we show here that the FRET dynamics 

can be modelled by Gaussian random processes with colored noise. Us- 
ing Monte-Carlo numerical simulations, the photon intensity correlations 
from the FRET pairs are calculated, and are turned out to be very close 
to those observed in experiment. The proposed stochastic description of 
FRET is consistent with existing theories for microscopic dynamics of the 
biomolecule that carries the FRET coupled dye pairs. 

Keywords: flourescence resonance energy transfer, single molecule, Ornstein- 
Uhlenbeck stochastic process, nanosecond dynamics, protein folding 

1 Introduction 

Recent significant advances in nano-technology make it possible to investigate 
molecular dynamics and structures at single-molecule level. Measurement of flu- 
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orescence resonant energy transfer (FRET) between a couple of dye molecules 
that are attached to complementary sites of a biomolecule like DNA or pro- 
tein is particularly useful, because the sharply distance-dependent dipole-dipole 
interaction between the dye pair can serve as a 'spectroscopic ruler' for the 
biomolecule ^ |21 El- FRET means a non-radiative quantum energy transfer 
from a donor that is a dye which initially absorbs light, to an acceptor which is 
another dye. 

The FRET can be considered in the framework of the theory of Forster pQ. 
An input laser light excites the donor, whose one of decay channels is to migrate 
its excitation energy to the acceptor via dipole-dipole interaction. The energy 
transfer typically finishes within nanoseconds. The requirements for FRET to 
occur efficiently are, at least, one of the chromophores should have a sufficient 
quantum yield and the donor fluorescence spectrum must overlap the acceptor 
absorption spectrum. 

Photon-photon correlations associated with the chromophores contain in- 
formation on the conformational distance of the biomolecule. Such informa- 
tion is usually washed out in traditional ensemble measurement, but is readily 
available in single- molecular measurement. Recently, using a Hanbury-Brown- 
Twiss 13] time-interval apparatus, Berglund et al. have measured photon 
intensity correlations for individual donor-acceptor pairs on DNA. To inter- 
pret the experimental data, they proposed a dual FRET model. A continuous 
model emphasizing overall conformational change of the biomolecule has also 
been studied by several authors jBl[7j. However, it is important to note that 
there exists non-trivial interplay between biomolecular diffusion, which under 
physiological conditions can change drastically the molecular conformation over 
nanoseconds, and the quantum optical processes of a pair of FRET coupled dyes, 
which are also of the time scale of nanoseconds. In this paper, we show that 
the laser-induced FRET dynamics can be modelled by the stochastic Ornstein- 
Uhlenbeck (OU) process [S]. The underlying stochastic process is a Gaussian 
random process [Hj with finite correlation time. The density matrix equations 
acquire the character of stochastic differential equations which can be solved 
using well-established methods. We shall demonstrate that the OU process is a 
good approximation to the FRET dynamics as measured on biomolecules. This 
paper is organized as follows. In section 2 we present our theoretical model 
of FRET. The idea of modelling FRET dynamics as an OU stochastic process 
is further elucidated in the context of biomolecular dynamics in solution. In 
section 3 and 4 we present the results of Monte-Carlo simulations and compare 
them with the experimental work. 

2 Theory 

In order to calculate the photon-photon correlation functions, we adopt a master 
equation approach. To derive the master equation for the system with two 
molecules coupled by FRET, we assume that the coherence time is much shorter 
than the time scale of experiments. We further note that a FRET coupled 
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system is not a cascaded quantum system j^, i.e., to obtain an equation for 
only donor variables, after tracing over the acceptor variables from the general 
master equation, is impossible. The master equation for the density operator ^, 
for the donor-acceptor pair couples four possible states \q)i (8) \p)2 here g = 0, 1 
{p = 0, 1) stands for the ground and excited states of the donor (acceptor) 
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where TrmiTi =1 — 5 are rate coefficients for the dominant processes: sponta- 
neous emissions of the donor and acceptor are described by the quantum jump 
operators Fi = cti (g) I2 and r2 = li ^ o'2; laser excitations of the donor and 
acceptor by T3 = I2 and r4 = ii ct^; and transfer of energy from the 

donor to the acceptor by Fs = (Ti ® 6-^. Here {a^) is the lowering (raising) 
operator between the excited and the ground states for the ith dye molecule 
{i = 1 for the donor and 2 for the acceptor). The master equation ^ ignores 
all coherence effects as they play no important role in FRET. The conditional 
photon count probability of the donor-acceptor system within a time delay r 
can be represented as a normally ordered correlation function |10 | 

(: Ut + T)i,{t) :) = Tr,,, {a+<7,f(r) {a,^(t)a+}} (2) 

here Q{t) is the stationary operator solution of Eq. and V{t) is the evolution 
operator of the whole system satisfying y(r)^(t) = Q{t + r). 

FRET measurements provide us information about dipole-dipole coupling, 
which varies as 1 /r^ ^I] . The FRET rate coefficient is given by 

TslO^rif^V (3) 



r{t) 

where i?o is the Forster radius, r{t) — \ri{t) — T2{t) \ is the distance between two 
chromophores. Considering intra-chain diffusion of a biomolecule, e.g., a protein 
molecule that carries the donor and the acceptor at a couple of complementary 
sites, the displacement of the donor-acceptor distance from its fixed (or initial) 
value, x{t) = rit) — r{Q), can be modelled as a Brownian motion in a harmonic 
potential. Since a biomolecule which consists of a large number of atoms and 
molecules, we essentially deal with the over-damped regime, where the Brownian 
motion is described by the Langevin equation [5] 

^ = -Xx{t) + \Ut) (4) 

where A = ijp /rj with lo being the frequency of harmonic oscillator, and rj the 
friction coefficient. The random force fx{i) is a pure Gaussian characterized by 
{fx{t)) = and {fx{t)fx{t + t')) = 2ed{t - t'). The solution of the Langevin 
equation, x{t), resembles a white noise: its equilibrium distribution is 
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and the correlation is 

{x{t)x{t')) =0Xe~^\*-*'K (6) 

Now turn to the FRET rate Eq. which can be rewritten as r5(<) = 
r5(0)/[l + a;(t)/r(0)]6 with r5(0) = Ti[Ra/r{0)]^ corresponding to the FRET 
rate at some fixed interdye distance. If the displacement x is small as compared 
to r(0), then we have 

ST.it) ^ x{t). (7) 

This is a possible interplay between the diffusion process and the FRET dy- 
namics, although neither they are completely uncorrelated nor fully correlated. 
It turns out that a variance (i.e., 'noise') of the FRET rate has the statistical 
signature of a white noise. A sign such as 'noise of the noise' usually leads to 
colored noise in the OU stochastic dynamics. Eq.(|7I) is a linearized approximate 
relation. Although, some dynamic details can not be seen for large variations, 
analytically, the exact numerical simulation may be done. This rather involved 
and so we have tried to produce reasonable results by retaining the leading 
term. This captures the physics reasonably well. Note that the OU stochastic 
process is stationary. Thus Eq. |^ becomes a stochastic differential equation as 
T^{t) = r5(0) +£,{t)- We assume that ^(t) is a colored noise, which is described 
by the Langevin equation |12| 

U^-K+\m (8) 



where time averages of white noise should be ri{t) = and ?]{t)ri{t') = 2D5{t—t'). 
As well known the Eq. ((S)), yields the steady state correlation function 

{mW)] = me-^l*-*'l (9) 

with ^(t) = and {...} denotes the stochastic average over the initial condi- 
tions ■ A parameter D might be proportional to the diffusion coefficient 9 
and A is the same in Eq.© and Eq.Q. The stochastic differential equation 
Eq.QJ will be solved using Monte-Carlo numerical simulations. A Box-Mueller 
algorithm and the Euler-Maruyama method have been used to realize the col- 
ored noise. Moreover, by virtue of an integral algorithm developed in we 
have verified that the Monte-Carlo generated correlation fits perfectly to its an- 
alytical expression Eq. (0. To achieve this, a stochastic averaging over as large 
as 1000 realizations is essential It must be borne in mind that the FRET 
rate is always positive which is done by keeping a background constant value 
r5(0). However, some large negative random numbers have to be omitted. To 
prove that these omitted random numbers do not play important role, we have 
also calculated the correlation for Tz{t) which decreases exponentially as given 
m Eq.Q. Namely, we assume that Th > r5(f) > where the FRET rate 
is finite. Vh (J^l) corresponds to minimum (maximum) inter-dye distance due 
to continues intrachain diffusion of the protein molecule. The quantities are 
determined by contour length and bending rigidity of the protein. In our case. 
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we have assumed that the FRET does not occur between distant dyes, so that 
Fi = 0. We take also Th ~ max{TQ{t)) as a maximum value of the generated 
random numbers. 



3 Results and discussions 

Using Monte-Carlo simulations we have calculated correlation functions as de- 
fined in Eq. ((SJ. In Fig. 1 we present the results for the normalized correlations 



Since only two parameters, D and A, are needed to completely specify an OU 
process, their determination would be a desired contact between theory and 
experiment. As we see in Fig. 1, the normalized correlation functions are very 
similar to those observed in the experiment For the data shown in Fig. 1 
the correlation time is taken to be Tc = 1/X = 7. In this case we take the 
parameter r5(0) to be 0.65 in order to ensure that the average FRET coeffi- 
cient to be around 1. The FRET coefficients fluctuate between = and 
Th = 5. The estimated Forster radius, for instance, for the TMR-Cy5 dye 
pair that is frequently used in biomolecular measurement, is about 53 A [2]. 
Given the calculated average FRET coefficient of 1, the average distance would 
be approximately 50 A. Intensity autocorrelations show the typical quantum 
feature of antibunching that is characteristic of emissions from individual dye 
molecules. Following the initial photon antibunching, photon bunching appears 
in the acceptor autocorrelation function. A sufficient stochastic deviation from 
its equilibrium distribution of the FRET rate is the hallmark of the generation 
of photon bunching, which is a tendency for clustered emissions. For large r the 
normalized correlation functions go to unity indicating uncorrelated emissions. 
The appearance of the photon bunching in the donor autocorrelation was first 
predicted theoretically by Haas and Steinberg A pronounced antibunching 
associated with a photon blockade effect and a photon bunching in the cross talk 
of the acceptor-donor pair, for short time, has also been discussed in [S]. We 
also notice in the experimental results that over longer times, photons emitted 
by the donor are becoming correlated with photons emitted by the acceptor 
and vice versa. It is also worth noting that because of off-resonant excitation 
of the acceptor, the corresponding rate for the acceptor is as small as /Fa with 
/ = 0.1. Otherwise, the acceptor would have already been excited by the laser, 
and FRET may not occur. For the same reason the laser excitation rate should 
not be too large. The bunching signal can be described approximately by an 
exponential y — 1 + Ce~^°'^ , with the fitting parameter Aq proportional to A. 
The tail of exponential decay, especially in the acceptor autocorrelation func- 
tion is determined by the correlation time Tc of the OU process. It is obvious 
that a larger value of the correlation time results in a longer tail. While being 
supported by the simulation results (Fig. 1), this point can be made clearer by 
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assuming that the intra-chain diffusion of the biomolecule that carries the dyes, 
and thereby fluctuation of the donor-acceptor distance is much slower than the 
quantum optical process of the donor and acceptor. The population on the 
excited states then adiabatically follows the slowly varying r5(f), so that the 
intensity of the acceptor can be approximated asj?] 

According to this adiabatic approximation in longer time delay r, the intensity 
correlation is an exponential [7] 

where we have assumed that the dual FRET rates result in a high value and a 
low value of the emission intensity, Ih and II , respectively. Using Th — 5 and 
Ti = obtained in the calculation for Fig. 1, we find that Ih — 0.93 and II = 
0.1 from Ea. (|ll() . The correlation function under the adiabatic approximation 
given by Ea. H12|) has been also plotted, see the dotted curve in Fig. 1. 



4 Conclusions 

We have shown how the FRET process as measured for biomolecules in so- 
lution can be modelled using the Ornstein-Uhlenbeck stochastic theory. This 
theory predicts the fluorescence intensity correlations from the FRET coupled 
dye pair, which are very similar to those observed in recent experiments. An 
analytic study based on the local linearization procedure also shows that it 
is consistent with existing theory for microscopic dynamics of the biomolecule 
that carries the FRET coupled dye pairs. The second-order intensity correlation 
functions for a FRET coupled dye pair, are largely determined by only a few 
statistical parameters of the FRET dynamics. It is found that the stochastic 
OU description helps elucidate the underlying mechanism for the experimen- 
tally observed fluorescence correlations that typically exhibit exponential decay 
over nanosecond timescale. 
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Figure 1: Normalized photon- photon correlation functions corresponding to 
donor-donor, acceptor-acceptor and donor-acceptor emissions. The parameters 
of the model (HJ are chosen as Ti^s.a = l,r4 = 0.1,r5(0) = 0.65, (T^it) ~ 
I7 < T5{t) < 5). The noise parameters are taken as Tc = 1/A = 7 and D = 7. 
The number of numerical realizations and a time step are taken to be iV = 1000 
and = 0.01, respectively. 
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